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Abstract 

The transition from toroidal counter- to co- rotation in the core 
plasma has been observed at L to H transition in several tokamaks. 
Spontaneous reversal has also been observed in TCV beyond a thresh- 
old in the density. We develop a model based on the following phe- 
nomenology: (1) the increase of the gradient of the pressure trig- 
gers formation on a fast time scale of cells of convection (similar to 
Rayleigh-Benard (RB), but with a single sign of vorticity); (2) poloidal 
rotation is induced by the envelope of the peripheric velocity of the 
convection cells; via the baroclinic term the gradients of temperature 
and density sustain the poloidal rotation against the decay due to the 
parallel viscosity; (3) the fast increase of poloidal flow induces a high 
time derivative of the radial electric field; (4) the neoclassical polar- 
ization creates a series of parallel accelerations (kiks on each bounce) 
of the trapped ions, leading to an increase of the toroidal precession 
or to its reversal; the source of energy is the work done by the radial 
electric field. (5) the diffusion transfers on resistive scale the toroidal 
momentum from the trapped ions to the untrapped ones. The cor- 
related interactions are examined and the estimated time scales are 
found to be compatible with the observations. 

1 Introduction 

The spontaneous reversal of the direction of the toroidal rotation in tokamak 
has been observed as a change accompanying the L to H transition or a slow 
variation of a plasma parameter, like the current or the magnetic field [1], 
[2], [3], [5]. Spontaneous reversal has also been observed in TCV beyond 
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a threshold in the density [B]. A review of these experimental facts can be 
found in Refs. [7] and [S]. 

There are several properties of the state of the system that seem to be cor- 
related with the transition consisting of the flow reversal. There is a strong 
radial gradient of temperature which is sustained by external input of energy 
into plasma. In general fluid configurations placed under strong thermal 
stress [for example in the classical Rayleigh-Benard (RB) system] the purely 
conductive state is replaced by a convective regime through generation of 
cells of convection |9]. The transport of heat and of momentum and angular 
momentum (via the Reynolds stress) become much more efficient compared 
with the diffusive regime. The onset of the convection structure takes place 
on short time scale compared with conduction or flow time scales. The event 
leading to formation of a regular convection pattern starts with the genera- 
tion of isolated hot and cold streamers (also called "plumes" or "thermals") 
that advance into the fluid environment [10]. The head of the streamer is 
deformed through interaction with the static fluid leading to "mushroom" 
shapes also known from Rayleigh- Taylor experiments. The streamer breaks 
down leaving vortices of two signs. These vortices are merging eventually 
creating the large scale flow, a mechanism which is universally associated to 
relaxation ^llj , ^12j . Under strong temperature gradient the whole process is 
very fast. 

The presence and/or the possible generation of closed, quasi- coherent 
convective flows in the meridional cross section of the tokamak plasma is a 
well investigated subject. The particular case of the zonal flows is included 
as a subset in this wide class of organized flows (for reviews see [I3], [H], 
|15j). Various mechanisms have been identified to be able to transfer en- 
ergy from the drift wave turbulence to the fiow in a convective structure 
extending over a substantial part of the plasma cross section. Modulational 
instability, wave-kinetic dynamical transfers in the spectral space, multi-scale 
interactions leading to negative viscosity are usually invoked [16], [17], [T8] . 
[T9j. A unifying and universal principle is the inverse cascade of energy in 
the two-dimensional turbulence [11], [20], [21], [22], [23]. Vortical structures 
consisting of robust and coherent fiow have been identified in regimes gov- 
erned by the Hasegawa-Mima and Hasegawa-Wakatani equations [15]. Due 
to the essential role of the vectorial nonlinearity the space scales of the these 
structures is initially at the level of the ion Larmor radius. The processes 
or merging lead further to creation of large scale structure. Another mecha- 
nism that has been studied is the tilting instability of eddies with creation, 
via Reynolds stress, of fiows over large spatial regions [21|. The structure is 
usually assumed with alternate sense of rotation, for several reasons. One of 
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them is the conservation of the angular momentum integrated over the space 
occupied by the cells of convection. This conditions is however imposed to 
closed systems, and even for them is not quite stringent [25] . A different rea- 
son comes from assuming the flow as a perturbation, expanding in Fourier 
series and retaining the lowest order which by its nature is periodic and leads 
to vortices with alternating sense of rotation [H] . It also provides compatible 
flows at the common boundaries of the neighbor cells. The latter argument 
should not be seen as a strong restriction if the distance between the rolls of 
convection is sufficiently large. Finally we mention a condition related with 
the space-time invariances of the equations, which appears to underlay this 
choice of the pattern of flow. Taking as a propotype the Hasegawa-Mima 
equation one notes the invariance to reversal of the sign of the vorticity, 
which implicitely imposes the presence in the solution of vortices of both 
signs. However it has been shown that besides the vectorial nonlinearity 
the scalar or Korteweg-deVries nonlinearity, of the type ip {dip/dy) where ip 
is the streamfunction and y the poloidal coordinate, is likely to be present 
in the equation [20], [2Z]- Although it does not influence the smaller scales 
where the instabilities are developing, it controls the large scales due to the 
lower differential degree. This breaks the invariance and allows solutions 
consisting of a set of rolls with unique sense of rotation. We argue that rolls 
of convection having the same direction of rotating flow are not excluded 
by any characteristic of the differential equations governing the quasi-two- 
dimensional tokamak plasma and are even more efficient for the transport of 
angular momentum in the meridional section. This pattern of flow has been 
derived |28] and it arises in many similar situations. In the science of atmo- 
sphere it is known that a distribution of azimuthal velocity in a vortex can 
be affected by perturbations leading to secondary vortices that move around 
the core of the parent system The geometry of the flow takes the form 
of a polygonal eye-wall, which means that a set of smaller vortices are placed 
symmetrically, with their centers placed on a circle around the eye-wall axis 
|30j . They all have the same sign of vorticity. Other examples are accretion 
disks [31] and the Ranque-Hilsch vortex tube [32] . 

In tokamak there is a poloidal rotation at equilibrium, due to the neo- 
classical processes. The poloidal flow is sheared since part of the plasma 
cross section is in plateau - banana and part is in Pfirsch-Schluter regime 
|33j . This just constitutes a seed which may favor a choice of the direction 
of rotation of the rolls that are going to be created. 

A source of flow that appears to be less investigated is the baroclinic term 
arising in the momentum equation of the ion component when it is trans- 
formed by a rotational operator into the equation for toroidally projected 
vorticity. It is well suited to exhibit local velocity shear that may accelerate 
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plasma flow, leading for example to generation of streamers. As shown below, 
the streamers constitute a basic agent for building up larger scale flows. 

It is well known that the RB transition from purely conductive to a con- 
vective regime takes place at a threshold, where the buoyancy force over- 
comes the viscosity [9]. In the present case the threshold has a different 
nature. The drive force in the early stage is the baroclinic term, dependent 
on the main gradients of the density and temperature, which are supported 
by the external input. Convective cells are likely to be created, as shown 
by numerous investigations (analytical, numerical, observations). The pat- 
tern should obey the symmetries required by the basic neoclassical rotation 
which favors formation of rolls all rotating in the same sense, equivalently, 
with the same sign of toroidally-directed vorticity. There is no restriction 
upon this choice coming from the conservation of the angular momentum in 
the meridional plane since momentum is carried away to the boundary (the 
system is open). Now, the peripheric velocities of the chain of large scale 
vortical cells constitues an effective poloidal rotation. The magnetic pump- 
ing opposes the poloidal rotation via the parallel viscosity which acts on a 
time scale determined by the ion- ion collision frequency [M]. This leads us 
to formulate a threshold condition: the gradients of density and temperature 
should be sufficiently high such as the baroclinic term to sustain the rolls 
and their induced poloidal flow against the parallel viscosity damping. The 
detailed examination of this threshold condition requires more quantitative 
analysis and will not be of the present concern. 

We develop a model based on the following phenomenology: 

(1) the increase of the gradient of the pressure triggers formation of cells 
of convection. The process is similar to the transition from purely conductive 
to convective regime in Rayleigh-Benard fluid but the rolls of convection have 
a single sign of vorticity; the baroclinic term drives a fast onset of this pattern 
of flow; 

(2) the peripheric velocities of the convection cells (with unique sense of 
rotation) coalesce into an envelope which is equivalent to a poloidal rotation; 
in usual circumstances the poloidal rotation is damped by the parallel vis- 
cosity (magnetic pumping) but this is now overcomed by the Reynolds stress 
in the regular pattern which is directly sustained by the gradients of density 
and temperature. These are, in turn, sustained by the external input. 

(3) the fast increase of poloidal flow induces a high time derivative of the 
radial electric field Er. The radial momentum balance would allow to divide 
the effect of an increase of the vg term to the two other terms, containing the 
diamagnetic and respectively the toroidal velocities. However the reaction 
of the diamagnetic velocity is slower since it depends on the change of the 
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density profile and the change of the toroidal flow velocity is mediated by a 
neoclassical effect. 

(4) the neoclassical polarization induced by the variation of the radial 
electric field creates a series of parallel accelerations (kick on every bounce) 
of the trapped ions, leading to an increase of the toroidal precession or to its 
reversal; the source of energy is the work done by the fast transient growth 
of Er supported in turn by the pressure gradient through the convection. 

(5) the diffusion transfers on resistive scale the toroidal momentum from 
the trapped ions to the untrapped ones, inducing the reversal of direction. 

The point (4) roughly describes this situation: on half of the banana (say 
in the same direction as the magnetic field B) the trapped ion feels a certain 
radial electric field and when it returns on the other half of the banana orbit 
(opposite to B) it already is acted upon by a different radial electric field. 
The key is that the change of the radial electric field and the ion bounce on 
banana take place on comparable time scales. This is possible due to the 
fast onset of convection and can be the explanation of the fast transitions 
observed in experiments |35] . 

Our model is similar to an " inverse Ranque-Hilsch" effect. In the Ranque- 
Hilsch vortex (not fully understood at this time) a strong rotation of a gas 
leads to temperature separation |^ and it is supposed that coherent vortical 
structures that breaks the azimuthal symmetry of the flow are responsible for 
the creation of the gradient of temperature [32] • In our case we dispose of the 
strong temperature gradient across the plasma section and the onset of cells 
of convection (which break the azimuthal symmetry) would create the other 
element, i.e. the poloidal flow. The process is sustained by the heat and 
density input which is converted into mechanical stress (radial variation of 
poloidal rotation) by the cells of convection. It is transitory with a fast phase 
(when the reversal occurs) followed by a slower variation which keeps the time 
derivative of the radial electric field at the level required for balancing the 
parallel viscosity damping. 

We focus in the present work on the logical consistency of the argument, 
examining the contributing physical processes under certain simplifications. 
The overall result confirms the validity of the model and justifies further 
study. 

The paper is organised as follows: we explain the role of the baroclinic 
term in sustaining the perturbations consisting of hot streamers advancing 
along the minor radius. The estimations shows that the time scale of the 
convection is very fast. Separately we provide an argument for the poloidal 
flow of the density toward the place of origin of the streamer, preventing a lo- 
cal depression of density and suppression or breaking up of the hot streamer. 
This argument is based on another manifestation of the baroclinic effect. 
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which in this case reveals the possibhty of the instabihty called Chaplygin 
gas, or anomalous polytropic. The next section is devoted to demonstrate 
that the peripheric velocity of the rolls of convection leads to an effective 
poloidal flow. The argument is here geometric and on this basis one can eas- 
ily superpose physical processes like Reynolds stress and coUisional diffusion 
of the momentum in the vecinity of the rolls of convection. The last part 
consists of simply adapting the arguments of Hinton and Robertson [37j to 
show that the neoclassical polarization induces an acceleration of the toroidal 
motion of the bananas. 



2 The hot streamers in the early phase of the 
onset of the convection rolls 

We consider a perturbation in the two-dimensional section of a tokamak 
plasma, consisting of the radial displacement of an element of plasma from 
the radius tq to higher radii, i.e. to toward the edge. For fast motion the 
element preserves its temperature and advances as a stream of hot plasma 
to higher r. The Figure 1 shows the geometry of the perturbation fixing also 
the system of reference [x = r,y = r6, z). The smooth surface represents 
the equilibrium pressure profile on which it is superposed the perturbation 
related to the advancing hot streamer. The curved arrow shows the effect of 
the baroclinic term. The two arrows at the right of the perturbations have 
the following meaning: the one directed against the perturbation represents 
the VT of the hot stream relative to the environment. The one poining 
approximately in the inverse x = r direction is the gradient of the equilibrium 
density (or mass), Vp. 

The density at the radius where the streamer begins its motion is n (ro) 
and the temperature at the same location is T (ro). The displacement of the 
hot plasma element along the radius is 6r = r — which is positive in our 
system. On this distance the equilibrium (background) pressure has varied 
and has the value 

dp 1 
p{r) = p (ro) + —6r = p (ro) - p (r) ri—rSr (1) 

where Lp < on an usual equilibrium pressure profile. This p (r) is the 
pressure in the environment plasma at a distance 6r from the point of initial- 
ization of the thermal streamer. Assuming that there is no loss of heat via 
diffusion the hot stream will have at the same position the pressure p(ro). 
The gradient due to this contrast between the hot streamer and the envi- 
ronment is directed (along poloidal direction at r) from the boundary of the 
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thermal streamer to its axis. We will choose to look at the left side of the 
hot streamer which means that the gradient of the pressure is directed along 
the positive y axis, eg. The width of the hot stream is of the order of the 
radius of the roll, R and the gradient of pressure transversal on the axis of 
the streamer is 

6p p (ro) -p{r) 11 

-^ = V B = VP{r)T^M- (2) 

Jri JrC l-'^pl 

In order to keep a realistic set of assumptions we must consider the loss 
of heat from the streamer via conduction to the environment plasma. An 
elementary option is to introduce an arbitrary nondimensional small factor rj 
and consider from now on the fraction 7] of the perturbation of the pressure. 
Now, the gradient of the perturbed pressure interacts with the equilibrium 
gradient of the density Pq to give the baroclinic term 

^Vpo xVp = \ |Vpol (-e,.) X |Vp| (ee) (3) 



Po Po 



Pi 



iVpol |Vp| 



and using Eq.(|2]) we have 



1 ^ ^ 1 1 p (r) (5r 

Po i^) \Ln\ \Lp\ R 

The two quantities: the depth of penetration of the streamer 6r and the 
radius of the cell of convection R produced by the roll-up of the streamer 
have similar magnitudes and we approximate 

Sr 

(5) 

The baroclinic term is directed along the negative z (toroidal) axis and en- 
hances the local gradients of the velocity of the streamer: the velocity on 
the axis of the hot streamer is accelerated. With this specification of the 
direction of the enhancement of the velocity by the time-derivative of the 
local vorticity, we will only take from now on its absolute value 

du 1 V 1 1 p(r) 

n 6tR min{r) ' \Lp\ ^ ' 

Here v = v,, is the speed of the hot streamer, along its axis. The interval of 
time bt will later be taken a characteristic turn-over of the convection roll. 
The gradient length of the equilibrium pressure is 

1 _\dp _ 1 (rj^dn ^^dT\_ 1 ^ndT 
Lp p dr nT \ dr dr J Ln p dr 
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with both terms having the same sign, negative along the increasing minor 
radius direction. From Eqs.dH]) and ([7]) we have the estimation 

I V T f I 



StR rrii \\Ln\ \Ln\ \Lt\ 

Now, if we assume that the gradient length of the density is much larger than 
{LtI, we will retain only the second term in the paranthesis 

— -'tn^lVTI (9) 

Taking 

5t ~ ^ (10) 

we obtain 

~ R^r] — -— |VT| (11) 



1 


1 


rrii 


\Lin 


1 ~ 


R 




V 




1 



rrii \L. 



n I 



The gradient of the equilibrium temperature at the position and on the spa- 
tial extension of the convection roll must be a fraction of the gross estimation 
T/R and we introduce an arbitrary nondimansional small parameter A 

|VT|~a| (12) 

it results 

- R'^^zAr-i^ = ^^^-M. (13) 

Except for the two small parameters rj and A this formula shows that the 
speed of rotation at the periphery of the roll is very high, comparable with 
the ion thermal velocity. Including a rough estimation for the conductive 
loss (similar to the Rayleigh-Benard case) by the factor rj = 0.1 and the 
variation of the equilibrium temperature over the cell's extension, by the 
factor A = 0.1, we still get a high value t> ~ 0.1 x vth,i- 



3 The local perturbation of the density along 
the poloidal direction 

3.1 The instability of the Chaplygin-type or negative 
pressure 

The perturbation is initiated at a radius tq, which means that the density in 
a small region around the circle tq will flow toward the point (say 9o) where 
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the streamer originates. There the poloidal aflux is converted into radial flow 
of the hot streamer. This radial flow is accelerated by the baroclinic effect 
and from the Eg. (1131) we see that the speed can reach high values. A local 
depression of density can occur if the rate of fllling with density the foot of 
the streamer takes place at a rate which is slower than that with which it 
is removed by the streamer. This is an important aspect which we examine 
in the following. The continuity equation with only the velocity along the 
streamer is 

dn d , ^ /X 

or 

Q Qij Q 

— lnn+ — ^ + fr7^1nn = (15) 
at or or 

The radial velocity Vr is zero for radii smaller than ro and rises to an almost 

constant value on a small radial distance ^ ~ Vr5 (r — ro) . The equation 

of continuity becomes 

— Inn = — (r — ro) (16) 
where ^ = -§i + Vr^- The solution is 

In n [r (t)] — Inn [r (to)] ~ / dt'vj.5 \r (t) — Tq] (l"^) 

where the integral along the trajectory is the inverse of the Lagrangian deriva- 
tion d/dt. The expression shows that the density is simply carried along the 
streamer, with a discontinuity in the small region at the foot of the streamer. 
This leads to a depression of density if the neighbor regions are not provid- 
ing aflux of particles. We will show that there is an instability consisting 
of strong modulation of the density on the poloidal circle ro which can lead 
to accumulation of density in periodic poloidal positions. Although this in- 
stability and the formation of the streamer are independent, they can act 
concurrently, avoiding suppression of the streamer or its breaking up with 
further propagation as blob. This is the reason for which we examine now the 
poloidal modulation of the density when a conjunctural factor (the streamer) 
imposes a flow along 6. 

Assuming incompressibility V • v = in the small region around the 
streamer's foot we have 6vg/6l0 = Sv^/Slr with comparable variation scales 
6lg ~ 6lr- The two velocities, vg and v = Vr are therefore related vg ^ v. 
Accordingly, the acceleration of the streamer velocity v along the radius 
induces an accelaration of the flow vg on the poloidal circle. 

^ ~ ^ (18) 
dt dt ^ ' 
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which simply means that the time variation of the velocity in the streamer 
is also the time variation of the poloidal dynamics. 

From now on our considerations will be refering to the poloidal direction 
only, more precisely to a narrow region around the circle of radius tq where 
the hot streamer is born. As said before the acceleration on the radial direc- 
tion of the flow in the hot streamer induces an acceleration of the aflux of 
particles, along the poloidal circle, dvg/dt, toward the foot of the streamer. 
In this region the dynamics along the hot streamer, essentially radial, is cou- 
pled to the dynamics along the poloidal circle. If the main characteristic of 
the hot streamer is the contrast of temperature relative to the environment in 
which it advances, the dynamics along the poloidal region regards essentially 
the density. This is because the temperature on the circle is constant while 
the particles must flow toward the streamer 's foot to sustain the streamer 
flux and this generates density perturbation. We introduce the density p 
on the circle and we have that in the region of perturbation induced by the 
streamer, p = Pq. The problem becomes one-dimensional, along the circle Tq, 
with the coordinate y = ro9 and we take the origin of y at the axis of the 
streamer, i.e. y = at ^o- For uniformity, we change the notation, Vy = vg. 

A new manifestation of the baroclinic effect takes place along the poloidal 
direction. We first remind the case of the streamer, where the baroclinic force 
arises from the contrast of temperature between the hot streamer and the en- 
vironment, in combination with the gradient of the background density. Now 
the baroclinic effect arises from the perturbation of density along the poloidal 
direction, in combination with the gradient of the equilibrium pressure, Vpo, 
in which the gradient of the equilibrium temperature Tq (r) is active. Assume 
that due to the flux of particles towards the foot of the streamer y — there 
is a perturbation in the form of a lump of density at this point,. We have, 
focusing our considerations on the left side of the hot streamer, dp/dy > 0, 



The effect of this baroclinic force is to enhance the flow along y toward y = 0. 
Taking a distance on y of perturbation of p (y) of the scale as the radius of 
the rolls of convection, R, we have an estimation of the effect on velocity 
coming from the left to y = 0, the axis of the streamer. There is acceleration 




]_ / dp\ dpo 



^\dy) dr ^ 



(19) 



dt dt dy dr 

and the balance of forces along y is 




dvy du _ jy\ ^ dp dpo 



Idp d{To/m,) 
p dy dr 



(20) 




Idp d{To/m,) 



dr 



) 



(21) 
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We set 



a = R 



d{To/mi 
Or 



R 



dr 



(22) 



with dimensions of squared velocity. 

We will express the term in the right hand side as a derivation to ?/, such 
as to exhibit the structure of the force resulting from the gradient of a scalar 
pressure, which we will have to identify. This means to find a function F (y) 
such as 



dp dy 



We obtain 



dFjp) 
dp 



-a 



-ap 



(24) 



and this allows us to characterize the term arising from the baroclinic effect 
force produced by the y-gradient of a pressure 



l^dp 
^ p'^ dy 



d__ 

dy^ 



(25) 



where the equivalent pressure, as a scalar function whose gradient is the 
force, is identified as 

p = —ap (26) 

and is negative. The equations of motion for an element of plasma situated 
on the poloidal circle at radius tq consists of the deformation displacements 
in the poloidal direction due to the force generated from the gradient of the 
scalar pressure 



dVy 

dt 



dVy_ 

dy 



1 d _ 

pdy 



and the continuity equation 



dp d 







(27) 



(2^ 



The same problem has been treated by Ott [38] where the temperature is 
negative. The pressure in the present case is negative with linear dependence 
on the perturbation of the density. This is the typical m = 1 case of the 
general treatment of Trubnikov and Zhdanov [SU] . 

In the Appendix A we show, following Trubnikov and Zhdanov [39], how 
this system is solved. The solution has the following characteristic: the 
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density on the poloidal circle has regions of broad, quasi-uniform magnitude 
separated by very high (spiky) accumulations. In the actual physical process 
the strong accumulations coincide with the initialization of the hot streamers 
and effectively support the flow of density they represent. 

3.2 Comment on the theory of the anomalous poly- 
tropic gas 

We comment on the applicability of the results obtained from the model 
of anomalous polytropic gas to the physical problem of breaking of the az- 
imuthal symmetry of the density and flow in tokamak with generation of 
radial streamers and further of rolls of convection. 

We have two distinct processes: (1) a hot stream of plasma leaving the 
region at radius Tq and advancing in the colder environment at higher r, 
sustained by a baroclinic force. The baroclinic force is generated by the 
combined effect of the contrast of temperature between the streamer and 
the environment and the gradient of equilibrium density profile. It is a fast 
process and a density depletion around tq can appear. (2) an instability of 
the density profile along the poloidal direction at Tq, again sustained by a 
baroclinic effect. This is determined by the combination between the gradient 
of the perturbed density along 6 and the gradient of the equilibrium pressure. 
The effect is generation of high accumulation of density at periodic positions 
along the poloidal direction. 

These two processes are interacting: the tendency of accumulation of den- 
sity at certain location on the poloidal circle supports the flow in the stream- 
ers. The decrease of the density at tq by the streamer flow is compensated 
by the aflux of density associated to the (Chaplygin gas - type) instabihty of 
the density on 6. Observations and numerical simulations suggest that the 
streamers ( "plumes" or "thermals") of hot fluid inhibit the instability on 
the transversal direction by the high rate of radial flow which sometimes lead 
to breaking- up of streamer from its original region and separation as a blob. 

The system of equations describing the instability in the early phase can 
be reduced to an ID model of a gas with anomalous polytropic (the pressure 
is negative and p ~ —cp). The solution consists of regions of broad and 
smooth variation of the variables (p, v) and very localised, spiky regions of 
high concentration of density. 

The instability of Chaplygin gas is fast: the streamers (which appear in 
the model as concentrations of density) are formed on a time scale given by 
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the equation 

- ~ ^1/2 (29) 



which results from the eigenmodes determined in the hnearized system (the 
dispersion relation already shows this characteristic scales). Here A and r 
are typical wavelength and time scales of the Chaplygin instability on y. Let 
us take pg = "^i'^ ~ 10^^ kg/m^, i? ~ a/10 ~ 0.1 (m), 

dPO 1 „^ n .r,A f J 



, , = 3 X lOM — (30) 

dr (a/10) / 

for n = 10^9 (part/m^) and T = 2000 (eV). Then 

^ 5 X 10^ (^) (31) 

Now we estimate the distance between the centres of two neighbor rolls, 
taking the pattern consisting of m = 5 convective rolls with their centres at 
tq ~ a/3, for a minor radius a = 1 (m), A = 27rro/m = 0.4 (m). Then the 
time scale is 

r = ^ ~ 10-^ is) (32) 

which is very fast. Moreover, this time scale is compatible with the time 
scale determined independently from the flow of the hot streamer, since from 
Eq.(fTOl) we also obtain r ~ 10"^.. 10"^ (s). 



4 The poloidal velocity as the envelope of the 
rotations in the cells of convection (rolls) 

The peripheric velocity of the convection rolls is transferred to plasma at 
larger r both by Reynolds stress and by collisional drag. On the other hand 
even the geometry of the flow around closed rolls leads to a continuous (but 
wavy) streamlines which can represent a continuous poloidal flow. We want 
to show how to obtain the envelope representing the connected poloidal flow 
arising from separated rolls having rotations in the same direction. There- 
fore we will restrict the treatment to the pure geometric properties of the 
ideal flow. The Reynolds stress and collisions will enhance the efficiency of 
producing a directed poloidal flow. 

The treatment is based on complex formalism [ID] for the plane flow. We 
consider the fluid in the complex z plane with a set of vortices placed at 
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different positions. We need the complex potential of the flow in the z plane, 
w (z). The problem is not solved directly in the physical plane of the flow 
but in a plane ( connected with z hj a conformal transformation. The idea is 
to take in the (^-plane a structure which is standard (canonical) . It consists 
of (1) the circle with the center in ^ = and unit radius |C| = 1 and (2) a 
set of M circles with centres Sj and radii Qj placed inside the circle \(\ — 1. 

By the conformal map ( ^ z {() this corresponds to M + 1 fixed circular 
obstacles in the z plane of the physical flow. Our rolls of convection are 
patches of rotating fluid for which it is only specified the circulation of the 
fluid around {i.e. the integral of the tangential velocity on the closed contour) 
7^ , j — 1,M. The fluid moves in the z plane with the condition of keeping 
the given value 7^- of the circulation at every circle j. This does not uniquely 
specify the flow. One needs the condition at infinity, far from the regions 
where the circular vortices are given. 

After defining the geometry it is possible to construct the complex poten- 
tial in the ( plane. It is easily shown that it is 



^(C) = -^ln 



(-a 



(c-d^) 



2 

— In C + const (33) 
27r 



This function represents the fiow around a unit vortex located ina E C— plane 
and is constant on the circle bounding the large domain \(\ — 1. The complex 

potential gives a circulation of —1 around the circle = j) y ■ d\ — —1 

ICI=i 

where v is calculated from w according to the formulas of the Appendix B. 

To simplify, we take the geometry in the 2;-plane of the physical flow as 
consisting of three convection rolls represented by three disks with centers 
at Z-i — —d, zo — 0, z+i — -\-d, with equal radii, s. The circulation around 
each vortex roll is 7. 

In the plane ( one has a circle of radius \(\ = 1. This will be transformed 
by the conformal map ( ^ z in the circle with center at 2^0 = 0, which is one 
of the flxed physical vortices (or rolls of convection), the central one. Two 
other circles are inside this |C| = 1 domain, located at S^i = —S, 5+i = +6 
and having equal radii = q. Now the geometry of the ^ plane is deflned. One 
has to write the complex potential in the ( plane. 

To impose a particular circulation on the contour of the j*^ circular patch 
of vorticity (or obstacle) from the set of M in the 2;-plane, one defines in the 
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(■-plane the potential 



(C, «) = -f In ( , , ""(ffi.^A for J = 1, M (34) 

The function that has been introduced is 

0, (C) = + (35) 

is a Mobius map. Then Gj (C, a) is the complex potential of a flow in the 
plane ( generated by a vortex of unit strength +1 in the point Q = a and 
making a circulation —1 around the circle 5j of radius q. This potential has 
zero circulation around any of the other M — 1 circles in the C, plane and 
around the circle |C| = 1- By superposition — 'yjGj{(,(3) , one 

j=0,j=l,...,M 

obtains in the plane ( the complex potential for the canonical geometry: |C| = 
1, S-i, S+i. Since we assume that the convection rolls have equal circulations 
7j = 7 for j = 1, 2, 3, we obtain w {() = -l^G- (C, «) - IqGq (C, a) - 
7_,_G'+(C,a). This complex potential is in plane and from here it must 
be carried back onto the physical (2;) plane, by a conformal transformation 
which requires numerical calculation. The resulting flow have streamlines 
that show formation of poloidal rotation as we depart progressively from 
the line where the centers of the vortices are located [lOj. Since this is a 
gemetrical property it means that the poloidal flow occurs simultaneously 
with the onset of convection rolls, i.e. on a very fast time scale. 



5 The effect of the fast increase of the ra- 
dial electric field and the acceleration of 
the toroidal fiow of bananas 

The high density and temperature gradients (or both), supported by exter- 
nal sources, lead to formation of large scale cells of convection. Due to the 
background neoclassical rotation the rolls have all the same direction of ro- 
tation. The baroclinic term induces high peripheric velocities in the rolls 
and the envelope of these peripheric flows is a poloidal rotation. The rise 
of the poloidal velocity takes place on the same time scale as the onset of 
the convection rolls and this is governed by the baroclinic term acting on 
the initial hot stream. It is a fast process, similar to the onset of convection 
in the classical Rayleigh-Benard system. One important aspect should be 
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underlined: the flow in the convection rolls (hence the poloidal rotation) is 
sustained by the gradients of the equilibrium temperature and/or the equi- 
librium density, which in turn are sustained by external sources. Therefore 
the flow in the rolls of convection (implicitely the induced poloidal velocity) 
is sustained by external sources. This completely changes the usual picture 
on the poloidal rotation, the damping due to the magnetic pumping is now 
balanced and overcommed by the external sources. This also suggests the 
existence of a threshold, which will be examined in a separate work. 

The dynamics in the radial direction consists of charge separation and 
has weak effect on the force balance along this direction. We consider the 
approximate radial equilibrium of forces 

Er^-^^ + vrBe - veBr (36) 
\e\ n ar 

as valid over the interval of time increase of vq. In Appendix C we argue that 
the fast change of vq is mainly reflected in a change of the radial electric field 
Er, on the same time scale, leading to the approximative balance 

Since we have 

it results that the time change of the radial electric field will modify the 
motion of the bananas. 

We are in the regime described by Hinton and Robertson [37] where during 
the motion of the particle on the banana orbit there is effective change of 
the radial electric field which acts on the particle. In the Ref.[37] the effect 
is found to be manifested in two ways. First, there is a steady radial motion 
of the bananas. Second, in each bounce period, there is an acceleration of 
the motion in the toroidal direction. The effects are due to the fact that 
(schematically) the particle feels a certain radial electric field on half of its 
orbit and a different radial electric field when it returns along the second half 
of the orbit. 

The main assumption that has led us to this estimation is that 

— - -- (39) 
dt~RSt ^^^^ 

i.e. the baroclinic term acts all around the circular path of the rolled up flow 
in the cell and accelerates the rotation. The time scale is given by the time 
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of rotation in the cell. 

The time variation of the poloidal velocity is the time of acceleration 



and the time derivative of the radial electric field is 

^ ' -veBr) ~ -V^^Bt (42) 



dt dt 

At every bounce the banana gets a toroidal acceleration [37] 

(dE \ 1 
~gf] -^^hounce (43) 

where 

Vi,th 

du ~ -r]X-—BT— = -Vi^thV^^-rn y^^) 

\-l^n\ -De Vi^th -De \-l^n\ 

This kick of acceleration, leading to the change of the toroidal velocity of the 
trapped particle 6u is due to the variation of the radial electric field during 
the time interval the trapped ion makes a complete bounce. The number of 
bounces that a trapped particle makes in an interval of time At is 

" = = ^^6) 

Otbounce Q-tto/Vi^th 

During the time interval At the total change of the toroidal velocity of a 
trapped particle is 

Au = nSu = -Atv^^j^nX^ (47) 

l-t^nl -De 

This change is substantial. On a time interval of one millisecond At = lO^^s, 
with Ln ^ 5 m , Bt/Bq ~ 5, Vi^th = 10^ m/s and taking as before r] = 0.1, 
A = 0.1 the change of the toroidal velocity of a single trapped particle is of 
the order of the thermal velocity. When the initial toroidal flow is opposite to 
this acceleration there is reversal of rotation. The energy requested for this 
change is provided, in our model, by the equilibrium gradients of temperature 
and density. 
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6 Discussion and summary 



There are other mechanisms that can sustain or modify the spontaneous 
toroidal rotation. The drift wave turbulence creates a statistical ensemble of 
fluctuating velocities and the coupling of the radial and poloidal components 
can have a non-zero correlation (Reynolds stress) with spatial dependence. 
The divergence of this correlation is a driving force for the rotation. Due 
to the long correlation length in the toroidal direction, the Reynolds stress 
is efficient in the poloidal plane and in particular is able to induce zonal 
flows. If it overcomes the magnetic damping it can induce poloidal rotation 
and so generate a non- vanishing dEj./dt and the subsequent neoclassical 
polarization. The latter modifies the toroidal precession, i.e. the plasma 
toroidal rotation. A detailed analysis of this connection should be made, with 
focus on transitory processes, as observed in Alcator C-Mode yj. However 
few aspects can be discussed qualitatively. Since it is based on a statistical 
ensemble the Reynolds stress needs a strong poloidal-symmetry breaking 
mechanism otherwise the average is zero. Indeed the symmetry breaking 
can exist, analogous to the fiows arising in the Rayleigh-Benard experiments 
after higher order bifurcations (fiows are called "winds"). However their 
manifestation is slower that all the other time scales: conduction, convection. 
It is difficult to see how the onset of a poloidal fiow via only Reynolds stress 
can take place on time scales less than milliseconds, if the background consists 
of random drift waves. On the other hand it is well established that large 
scale convective fiows have much higher Reynolds stress than the 2D quasi- 
isotropic turbulence. This is partly the reason for the present investigation, 
where we precede any effect of change of the rotation by the onset of a 
chain of convection cells. We further note one aspect that results from the 
combined effect of convective cells and the Stringer drive and compare it with 
the Reynolds stress. A chain of convective cells (with either unidirectional 
or alternating vorticity direction) induces a periodic modulation of the local 
rate of transport. The neoclassical Stringer effect converts this modulation 
into an effective poloidal torque, possibly higher than the damping due to 
magnetic pumping. In a similar case but in the absence of convection cells 
|41j . the ratio of the magnitude of the Reynolds stress over the Stringer 
drive was found small. The presence of convective cells enhances both the 
Reynolds stress and the Stringer drive. Further study is necessary for a 
correct comparison of magnitudes. 

We have proposed a phenomenological model for the reversal of the 
toroidal rotation in tokamak. The model assembles effects that separately 
are known to appear in tokamak plasma and have been studied both exper- 
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imentally and theoretically. The generation of convection structures with 
flows that are quasi-coherent on large scales in the meridional section is in 
general admitted as a particularly effective way of transport of angular mo- 
mentum and of energy across the magnetic field. We consider the situation of 
cells of convection with the same sense of rotation flow, disposed around the 
magnetic axis in the region extending from the core plasma to the confining 
zone, a structure that is compatible with the basic neoclassical flow. It is 
also inspired by atmospheric and protoplanetary flows where such distribu- 
tion occurs in hydrodynamic regimes very similar to the one in tokamak. The 
convection is sustained by equilibrium gradients of density and temperature 
(which in turn are fed by external sources) and this drives the poloidal ve- 
locity which is generated by the combined, unidirectional, peripheric flows of 
the cells. Then the decay due to parallel viscosity (magnetic pumping) can 
be overcomed. The time scales are short due to the baroclinic term which 
accounts for the generation of vorticity. The fast increase of the poloidal ve- 
locity induces a fast variation of the radial electric fleld, on a time scale that 
is comparable with the bounce time. On part of the banana trajectory the 
ions see a magnitude of the radial electric field and on rest of the trajectory 
they see a different magnitude. This produces both radial motion and the 
change of velocity in toroidal direction. The collisions will spread this drive 
to the rest of the ion population. 

The principal merits of this theoretical explanation is that it provides 
a connection between poloidal and toroidal rotation. Usually the direction 
of influence is reversed and based of instabilities: the sheared parallel flow 
induces instabilities whose turbulent Reynolds stress acts to generate and 
sustain the poloidal rotation against the magnetic pumping decay |12], [13] . 
This mechanism is active in situation where the turbulence is known to induce 
poloidal flows acting as internal transport barriers, while the basic flow is 
observed to be toroidal. It can interfere with the generation of convective 
structures, for example in formation of Goertler vortices [44j, etc. In the 
present case the poloidal flow is a geometric consequence of the flow in the 
convection structures. The driving effect of the poloidal flow on the toroidal 
flow has a fast transient phase which is decissive since it involves neoclassical 
polarization on the time scale of the bounce on bananas. 

The different components of this phenomenology appear to be compatible 
and to combine in a structure which is logically coherent. Future work is 
needed for a more detailed examination of this theory. 

Acknowledgement 1 Work supported partially by the Contracts BS-2 and 
BS-I4 of the Association EURATOM - MEdC Romania. The views presented 
here do not necessarly represent those of the European Commission. 
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A Appendix. The Chaplygin gas instability 

A.l The formalism for the Chaplygin gas instability 

We will linearize the equations (127|) and (128|) . i.e. 

| + |W = (A.1) 

and 

dVy dVy Id / . X 



around the state p = , Vy = and obtain 
d ^ a 



7^(Po + P) + 7T-[(Po + P)^.] = (A.3) 



We have 



dvy dvy Id 

+ = ~^P^ 

dt ' dy Po + pdy 



dp dvy , . , , 



Po V Po/ Po 

Now we take p ~ exp [i {ky — ut)] and Vy ~ exp [z {ky — ut)] and inserting in 
the linearized equations we get 



— iujp + p^ik 



ikp 

lUJ pQ . 



(A.5) 



or 

a/2 



uj = ±ik{ay^" (A.6) 

The positive imaginary part indicates instability, similar to the Rayleigh- 
Taylor case. 

The numerical solution of the system of equations [38j shows that there 
are broad regions of weak variation and very peaked accumulations, which 
evolve to singular spatial spikes. 

The equations in Ref.[39] page 143 are relevant to the present case, for 
our negative pressure 

p = -ap (A.7) 
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If we introduce this expression in the equation of motion we obtain 



dVy dVy Id 1 9 , / . „S 

at ox poy Pq + poy 
and the system can be rewritten in terms of two variables p and Vy. 

dvy dvy a d ^ 
dt dy po 9y 

The equations belong to the elliptic type and therefore they do not have 
running waves. These equations correspond to the case m = 1 in the general 
treatment made by Trubnikov and Zhdanov and it is shown that they lead 
to the Laplace equation. 

We now try to apply the hodograph transformation to the equations of our 
case. Instead of the variables 'p{y,t) and Vy {y,t) we will use the inverted 
variables y (p, Vy) and t (p, Vy) with the relationships between the derivatives 

(A.IO) 



dVy 


]_dy 


9Vy ^ 


1 dt 


dt 


r\ — ■ 7 

w op 


dy 


w dfp 


dp 


1 dy 


dj) 


1 dt 


dt 


w dvy 


dy 


w dvy 



The Jacobian of the transformation is 

dy dt dy dt 

dp dvy dvy dp 

Now we substitute these equations in the system, to make the transfor- 
mation to the inverted variables {y,t) instead of (p,Vy). 

1 dy f I dt\ 1 dt 1 dt 



w dvy \ w dp J w dvy \ w dp 

We separate the derivative dy/dvy 

dv dt dt 

= -(P. + P)^ + «,5- (A.13) 



21 



Making the same substitution in the second equation we have 



dvy 
dt 



+ u 



a 



d 



w dp 



^ dy 
w dp 



-P 



a 1 di 

Po w dvy 



(A.14) 



and we separate the derivative dy/dp from this equation 

dy dt a dt 
dp ^dp pgdvy 



(A. 15) 



Now we have two partial derivatives of y to the two new variables (p, Vy). 
We impose the condition of compatibility of the mixed second order deriva- 
tives 

^ ^ (A.16) 



dp dvy dvy dp 



which leads to 



1 d 
pQ + 'pdp 



(po + py 



dt_ 
dp 



+ 



a_dH 
p^dvl 







(A. 17) 



To go further (and show that this equation leads to the Laplace equation) 
we make another substitution of variables 



r = {po + p) 



1/2 



2 a/ct/Po 



Then 



Vy — az 



where a = 2 \J (J j p^ and the change of variables implies 



d d dp ^ d 

dr dp Or dp 



d_ 
dz 



d dvy 
dvy dz 



a 



d_ 
dv 



or 



d _ 1 d 
dp 2r dr 

Inserting in the equation, we have 



1 



d 



Po + pdp 



(po + py 



_d_ 
dv„ 



df 

dp 



ld_ 
a dz 



^dH 

pQdvl 



(A.18) 



(A.19) 



(A.20) 



(A.21) 



(A.22) 
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which can be rewritten 



^ + ^^ + ^-0 (A 23) 

Following the general procedure (of Trubnikov and Zhdanov) we make 
the substitution 

^ = rt (A.24) 

and replace the variable t hj ip 

t = ^ (A.25) 
r 

dt 1 ^ Idip 

dH _ 2 ^ 1 di/j 1 d^j ^1 



and we obtain 



^ + -^-^2^ + ^ = ^ (A-26) 
This is the Laplace equation for a potential of the form 

(r, yj, 2;) = (r, 2;) cos (A. 27) 

A. 2 The solution of the Laplace equation 

In Ref . [39j it is explained the so-called " Evolutionary principle " , a restric- 
tion placed on the solutions of the Laplace equation: the physically useful 
solutions must satisfy the condition that, returning in time to the infinite 
past t — )• —00 the function \l/ or (p, v) must take the values of the equilib- 
rium profile. This means p ^ Pq and Vy — )■ 0, or, equivalently, r = 1 and 
z = 0. In accord with the last form of the Laplacean expressed in variables 
{r,(p,z) this corresponds to a circle or radius 1 in the plane z = 0. All 
the initial conditions at infinite past must be taken on a thin region around 
this circle. Since the Laplace equation also describes the electrostatic field, 
Trubnikov and Zhdanov exploit this analogy by formulating the problem as 
the solution of an electrostatic problem in space, generated by a certain dis- 
tribution of charges on the circle r = 1 and z = 0. The solutions that are 
admissible must be singular on this circle. The distribution of charges is 
equivalently described in terms of multipoles. Trubnikov and Zhdanov show 
that only two multipoles are useful for applications: the Coulomb distribu- 
tion and the dipole distribution of charges on the circle. The first provides 
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periodic solutions in real {i.e. physical) space. The " 2-multipole" (dipole) 
generates three types of solutions in the physical space: the hill, the well and 
the doublet consiting of well+hill. 

The solution of the Laplace equation in the space (r, ip, z) with boundary 
conditions given on the circle (r = 1, 2; = 0) requires to assume the toroidal 
geometry and therefore to introduce the system of coordines 

r = Asinh,^ , 2; = A sin 77 (A. 28) 

where A = (cosh^ + cosr/)"^ and the metric is 

cir^ + dz^ = {XdCi^ + {Xdrjf (A.29) 

The Lame coefficients are — X , — r — Asinh^ and — X. The 
Laplace equation is 

In our particular case, which corresponds to m = 1 of the general theory, 
we must seek solutions of the type = (r, z) cos ip and after the substitu- 
tion 

^{r,z)^r-^''^f (A.31) 
the amplitude is written as a series of terms with separated variables 

f-J2f-(0^os{nrj) (A.32) 

n=0 

Now it is introduced a new variable instead of ^ 

X = cothe (A.33) 

and this transforms the Laplace equation into the equation for Legendre 
polynomials 



d 

dx 



dfr, , . . n 



+ 



2 



fn = (A.34) 



for i/ = —3/2. The solutions are, in general 

P;(0 and Qli^} (A.35) 

However only the solutions (C) must be retained since they have the prop- 
erty to become singular on the circle r = latt— )■— 00. 

The solution is characterised by large regions where the function is smooth 
with weak variation and local spikes which evolve in time to singularity. The 
distribution of density is then compatible with the independent creation of 
streams of plasma which require strong density flow. 
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B Appendix. The complex potential of flow 
of a system of circular vorticity patches 

The flow is described by a complex function named complex potential w. It 
is composed from: 

1. the real part, (f{x,y), is the potential of the flow, i.e. the function 
from which the components of the velocity are derived 

This is the simple expression of the physical fact that the flow is poten- 
tial, i.e. it is derived as the gradient of the scalar function v = V</?. The 
flow is therefore irrotational, the vorticity is zero u? = Vx (V</?) = 0. 

2. the imaginary part, ip {x,y) which is the streamfunction and is related 
to the potential function ip by the fact that it also gives the components 
of the velocity 

dip dip 

Then 

d^^_dip_ d(f^di^ 
dx dy ^ dy dx 

which are the Ricmann-Cauchy conditions for holomorphy of the com- 
plex potential w (x, y). The lines of streamfunction ip (x, y) are perpen- 
dicular on the lines of the potential (p. The velocity is tangent to the 
lines %p (x, y) =const. 

Then the complex potential is 

w{x,y) = (f {x,y) + iip {x,y) or (B.4) 
w {z) = ip {z) + lip {z) 

Since (p and %p are the real and imaginary part of a holomorphic function (in 
a region of holomorphicity) both are harmonic functions 

A(/j = (B.5) 
Az/) = in the regions outside vortices 

and the flow has no vorticity with the exception of the closed regions where 
the circulations is defined, 7^ . 
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C Appendix. The fast increase of the poloidal 
velocity 

The problem consists of calculation of the radial polarization current when 
the poloidal velocity is strongly forced. The process has two stages: the first is 
the onset of the poloidal rotation, as envelope of the peripheric flows of a chain 
of convection cells in the meridional plane; the second is the sustainement 
of the poloidal rotation against magnetic damping and viscosity. The first 
stage is very fast, as is suggested by the similarity with the Rayleigh-Benard 
bifurcation conduction/convection and as is confirmed by the evaluation of 
the drive of the streamers by the baroclinic term in the vorticity equation. 
In this stage the rise of the poloidal flow generates a time-varying radial 
polarization electric field which acts on trapped ions on the time scale of 
their bounce producing a toroidal acceleration. The work done by the radial 
electric field on trapped ions is supported by the energy coming from the 
flow in the convection cells and this one is in turn supported by gradients of 
equilibrium plasma parameters (ultimetely from external input). 

In general the radial projection of the ion momentum connects terms that 
contain the poloidal velocity f*, the toroidal velocity f*, the radial electric 
field Er and the diamagnetic term ~ dpio/dr. Since now the poloidal velocity 
is forced, it is important to see which of the other terms will get the maximum 
effect on the short time scale of the onset of the convection rolls. Based on 
reasonable assumptions we show that the drive of dVy/dt goes primarly to 
the radial electric field {dEr/dt, the polarization). This further justifies to 
look to its effect on trapped ions. 

The estimations of the important effects will proceed along (and it is 
part of) the derivation of the initial value problem for the poloidal velocity, 
as shown in Ref. The radial current results from the Ampere's law 

dE 

V X B|^ = = fiojx + ^of^o-Qf (C.l) 

where 

J, (t) = -envl + envl (C.2) 

We use a slab geometry where the coordinates are: x = r with e^. = 
—Gr (toward the plasma interior), y = rO with = (poloidal), and 
z = Rip with Gz = (toroidal). The conservation of the ion momentum 
in the radial direction is 

'dt ^ '"y^J ^ ~~Bt ^ ^^'^'^ ^ en.vlB, - eriivlBy + (C.3) 
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The lines of flow in the poloidal direction are wiggly due to the discrete chain 
of rolls but on the average the term with the poloidal {y) derivative of the 
radial velocity vanishes. In the following we we neglect the resistive terms. 

^ eBz dt Bz Bz ^ UimiQi dx 

the last term being the diamagnetic flow contribution. The equation for the 
poloidal velocity f * is indicative for the effect that can result from the fast 
time variation, when the poloidal flow is strongly forced by convection rolls. 
The time variation of Vy can induce a time variation of the diamagnetic term, 
but this can be considered small since the change of the equilibrium proflle 
of the density involves diffusive time scale. The change of the density due to 
the perturbation {i.e. the convection cells) appears here in second order and 
will be neglected. 

Further, we have to bear in mind our main objective, which is to show the 
possibility of coupling of the poloidal and toroidal rotations via the neoclassi- 
cal polarization, acting on trapped ions. For this reason we can postpone the 
inclusion of the direct {i.e. not mediated by neoclassical polarization) effect 
of dVy/dt upon the term ^vl. The relative importance of this term will be 
discussed at the end of this Appendix. Then we will use the approximations 



The conservation of the ion momentum in the poloidal direction is 

'dt ^ ) " ^ ^""'^^ ~ '^'"'''^^^ ^ ^^-^^ 

Since this calculation is made on the points along the poloidal circle where 
the convection rolls are closed the radial ion velocity v^. on this circle and 

the nonhnear term are small. Also the poloidal {y) components of the 

electric field Ey and of the gradient of the pressure are vanishingly small 
compared with the rest of the terms (actually they are neoclassical correc- 
tions related to trapping and toroidal viscosity, which are here comparatively 
small). It results 



V dt 

Now we use the Ampere's law 



1 + 7^ 



1 / dv. 



ili \ dx 



(C.7) 



dEx 

euevl = -euivl - £o-^ (C.8) 
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In this expression we insert the radial electric field from Eq. flC.5p and the 
radial ion's velocity from Eq.f lC.7p . We obtain a single equation for the 
poloidal ion velocity 



en 



eoTTli 



efli \dt^ J U eQf \ dt^ J \dxdt J 
enf \dt J \dxdt) f/3 efif \ dt J \dxdt 



1 

LP 



1 / dv^ \ 

with the notation U = 1 + ^ ( -^^ j . For the phase of the acceleration the 
variation on x (radial) is not relevant and a simplified form is obtained 

The calculations must be repeated for electrons, to derive the expression of 
The relevant terms have however factors of and since Qe ^ the 



V, 



electron current can be neglected 

The quantity e± = eo{l + c^/v\) is the effective plasma dielectric in the di- 
rection transversal to the strong magnetic field. The initial condition is the 
acceleration of the poloidal fiow resulting from forcing at the onset of the 
convection, K = dVy/dt. As said before, this process (similar to RB bifur- 
cation) is expected to be very fast but its details are still beyond analytical 
description. 



where G (t) is the Heaviside function. With the time scaled t tQc the 
equation becomes 

w + bw = -KB (t) 

Here w (t) = dVy/dt and the dot means derivative to the scaled time variable. 
The notations are h = e±/eo and K = fl^^K. The solution, obtained by 
Laplace transform [15] is 



-K ( -K 



w{t) = — Y — ^(""^0 T~ ) cos {i^i/^o)^^'^ t { — — ) sin (e^/^o)^ 



1 \ 



1/2 



— -, — \-\ Wo 1— COS [e±/eo) t +wi — — 

28 



For simplicity we have not included the terms that induce the decay of the 
oscillations [H]. The initial conditions are given for the first and the second 
derivatives of the velocity, respectively wq = w (t = 0) , and wi = w (t = 0). 
From the definition of K we have Wq = w {t = 0) = ^~^K = ^ (j^^ 

The non-oscillatory part of acceleration w (t) = results from estimations 
made in the Section 2. During the onset of convection cells the poloidal 
velocity reach values comparable to the ion thermal velocity ~ 10^ (m/s) in 
the characteristic time of advancement of the streamer, 10^^ (s). Then the 
initialization is {dVy/dt)^^^ ~ 10^^ {m/s"^) and the non-oscillatory part has 
the order of magnitude 



dVy 



dt J e^/ £q 



10« (m/s^) (C.12) 



Now we can return to Eq. (lC.4p and use this estimation to compare the terms 
that are supposed to share the effect of the fast increase of the poloidal 
velocity. The only term that may preoccupate contains the toroidal velocity. 
Taking as reference the experimental observation [2] the toroidal velocity is 
of the order of 10^ [m/s) and any change takes place on a time scale of 
fractions of a second. This means that the observed toroidal acceleration 
is substantially slower than Eq. (lC.12|) . Besides, this term is multiplied by 



the small ratio By/Bz ^ 1. We conclude that the fast time variation of 
the poloidal velocity can only be followed by the radial electric field, which 
justifies Eq. (lH7ll . 
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Figure 1. 

Figure caption: The system of reference and the geometry of the per- 
turbation. The smooth surface represents the equilibrium pressure profile on 
which the temperature perturbation of the hot streamer is superposed. The 
curved arrow indicates the direction of the baroclinic velocity acceleration. 
The two arrows at the right side of the perturbation represents the gradients: 
the arrow directed to the left is VT of the hot streamer; the arrow directed 
to the negative x axis is Vrio of the equilibrium density. 
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